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ABSTRACT 

We discuss a pulsar acceleration mechanism based on asymmetric neutrino emission from the direct quark Urea process in the interior 
of proto neutron stars. The anisotropy is caused by a strong magnetic field which polarises the spin of the electrons opposite to the field 
direction. Due to parity violation the neutrinos and anti-neutrinos leave the star in one direction accelerating the pulsar. We calculate 
for varying quark chemical potentials the kick velocity in dependence of the quark phase temperature and its radius. Ignoring neutrino 
quark scattering we find that within a quark phase radius of 10 km and temperatures larger than 5 MeV kick velocities of 1000km 
s" 1 can be reached very easily. On the other hand taking into account the small neutrino mean free paths it seems impossible to reach 
velocities higher than 100km s -1 even when including effects from colour superconductivity where the neutrino quark interactions are 
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1. Introduction 



Ten years after its disco very in 1967 Wyckoff and Murray 
dWvckoff & Murra"vlll977l) estimated for the Crab pulsar a tan- 
gential velocity of 123 km s _1 out of the centre of the supernova 
remnant. This high transverse speed is not very different from 
the value of 140km s considered today dHobbs et al.l 12005), 
nevertheless compared to other pulsars the Crab pulsar speed is 
still quite low. Pro per motion measurements of 233 pulsars by 
iHobbs et al.l d2005l) imply a transverse velo city spectrum which 
is rang i ng up to > 1000km s " 1 (see also ICordes & Chernofn 
dl998l) ). IChatteriee et alJd2005l) found a directly measured speed 
for a pulsar of v = 1083^g 3 km s _1 . They furthermore argued that 
a natal kick is required to impart the measured proper motion. 
The reconstruction of the initial veloc ity distribution of these fast 
pulsars is still an issue of discussion. lArzoumanian et al.l d2002l) 
found a two-component velocity distributio n with characteris- 
tic vel ocities of 90km s _1 and 500km s _1 . iBombaci & Popovl 
d2004l) proposed that the simultaneous existence of neutron stars 
and quark stars could explain the bimodal distribution. On the 
other hand, IHobbs et al.1 d2005l) find that the velocities of young 
pulsars are well fit by a single Maxwellian distribution with a 
mean velocity of 400km s _1 . Another interesting point concern- 
ing pulsar velocities is t heir direction. By fitting pulsar wind 
tori iNg & Romanil d2004l) obtain nearly model independent es- 
timates for the neutron star spin orientations which can be com- 
pared with the axes of proper motion. For the Crab and the 
Vela pulsar a good alignme nt between the veloc ity vector and 
the rotational axis is seen dBrisken et alj 120051) . In addition, 
iJohnston et alj d2005l) examined the position angle of the rota- 
tional axis due to polarisation measurements of 25 pulsars and 
concluded that there is a strong indication for an alignment of 
the velocity and the spin vector. 

Since the discovery of the high neutron star space veloci- 
ties there has been a large number of suggestions to explain 
their origin. The general assumption is that during a certain 
time period in their evolution neutron stars must have experi- 



enced an accelerating kick, a so-called pulsar kick. An overview 
of the most common kick mech anisms can be found in e.g. 
dWang et al.ll2006h lLai et al.ll2001l). The hydrod ynamical mech- 
anism dBurrows et al.1 120061 : iJanka et akl 120071) is based on an 
asymmetric supernova explosion where the explosion is stronger 
in one direction and the neutron star receives a recoil. Kick ve- 
locities of more than 1000km s" 1 are p o ssible for s ucces sful 
supernova explosions dScheck et al.ll2006l) . lGott et alj d 19701) on 
the other hand considered the disruption of a binary system due 
to the explosion of one partner where the remaining star leaves 
with nearly its orbital velocity. Another possible acceleration 
source can be asymmetric low-frequency electromagnetic radia- 
tion due to a n an off-centred rotating dipol e (electromagnetic ally 
driven kicks dHarrison & Tademarul 1975b ). The combination of 
a rapidly rotating collapsing iron core and a magnetic field can 
lead to the formation of two jets at the poles. In the case of an 
asymmetry betw een these jets the supe rnova explosion can lead 
to a large kick dKhokhlovet al.lll999l) . Neutron star accelera- 
tion can also be based on the fragmentation of a rapidly rotat- 
ing core into a double proto-neutron star, where the explosion of 
the lighter proto-neutron st ar could accelerate the remaining one 
dColpi & Wassermanll2002l) . 

A widely discussed class between acceleration mechanisms 
for pulsars is base d on asymmetric neutrino emission 
dDorofeev et alj 119851) . Neutrinos can be a source for pulsar 
kicks due to the large energy release during the supernova and 
the following proto-neutron star evolution. A neutron star mov- 
ing with a velocity of 1000 km s _1 has a kinetic energy of ~ 10 49 
erg while the energy released in neutrinos is much larger, around 
~10 53 erg. From momentum conservation one can see that an 
asymmetry of about 3% could accelerate a neutron star of 1.4 
M Q to velocities of 1000 km s _1 . The momentum of neutrinos 
(p v ) and the neutron star (p ns ) reads: 



Pv = 



c 3 ■ 10 10 cm/s 

Pns = M ns ■ Vkick 



3 ■ lO^erg _ 1Q43 erg-s 



cm 



(1) 
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= i.4m ■ 10 8 — - 2.8 ■ 10 41 ^^ - 0.03 Pv 
s cm 



(2) 



Considering the violent birth of neutron stars as well as their 
large magnetic fields one could assume that anisotropics in 
the neutrino flux sho uld be rather standard than an exception. 
ISocrates et al.l d2005l) discuss a radiatively driven magnetic in- 
stability which leads to a local neutrino lu minosity enhance- 
ment during the Kelvin-Helmholtz cooling. IHorowitz! d2005l) 
discusses neutrino flux anisotropics stemming from parity vi- 
olatio n in weak interactions in prese nce of a strong magnetic 
field. iHorowitz & Piekarewicz (1998) examined amongst oth- 
ers anisotropics in neutrino-momentum distribution from cap- 
ture of polarised electrons, whereas the resulting asymmetry 
of the neutrinos originates from the difference of the nucleon 
weak axial and vector coupling constants. The required mag- 
netic fields were found to be in the range of ~ 10 16 Gauss 
throughout the whole proto neutro n star. In case of smaller mag- 
netic fields IHorowitz & Lil (fl998) studied cumulative parity vi- 
olation in neutrino elastic scattering from polari sed neutrons in 
proto- neutron stars. However, it was shown by iKusenko et alj 
(119981) that for such processes no asymmetry would be gener- 
ated in thermal equilibrium even in the presence of parity viola- 
tion with a nisotropic scattering ampl itudes (kn own as the "no -go 
theorem'"). lKusenko & Segre]Il996l) as well as iGovalld 19991) ar- 
gued that the tau neutrino sphere is located at a smaller radius 
than the one for electron neutrinos, so that electron neutrinos 
can oscillate to r-neutrinos in between the two spheres due to 
the Mikheyev - Smirnov - Wolfenstein effect (MSW-effect) and 
escape the proto-neutron star. The distortion of the conversion 
surface due to a strong magnetic field would lead to tau neu- 
trino emission from different temperatures at different regions 
and therefore to a pulsar kick. A method to avoid the prob- 
lem of small neutrino interaction rates completely involves ster- 
ile neutrinos. In contrast to normal neutrinos whose anisotropy 
is washed out for high interaction rates sterile neutrinos would 
keep their asymmetric distribution due to the vanishing interac- 
tion rates. Sterile neutrinos can be produced due to the MSW- 
effect from normal neutrinos or directly in we ak processes sup - 
pressed by the square of the mixing angle (see lKusenkol d2005l) : 



brid stars form a new stab le solution of the Tolman- 
Oppenheimer-Volkoff equation dGlendenning & Kettner i [200C : 
iMishustin etaT] 120031: ISchertler et alj 120001: iFraga et alj I2OOI 



Fryer & Kusenko (2006) and references therein). 



Ng & Romanil (120071) recently explored the correlation between 



pulsar spin and proper motion assuming a single acceleration 
kick during the proto neutron star cooling and find that the thrust 
scaled proportional to the neutrino luminosity. They find that the 
preferred fit parameters are consistent with a magnetic field in- 
duced asymmetry of a neutrino driven kick where the character- 
istic timescales for the anisot ropic emission are 1 -3 s. Similar 
conclusions were derived by IWang et ail d2006l) who found a 
kick timescale of hundreds of milliseconds to 1 s which fits well 
with magnetic-neutrino driven acceleration. 
In this work we investigate such a neutrino kick mechanism for 
proto-neutron stars with an exotic quark matter core. Besides 
the traditional neutron star there are various predictions on the 
presence of exotic matter in the neutron star interior, as hyper- 
ons, Kaon conden sation or quark matter (for an overview see 
e.g. IWeberl d2005l) ). which is expected to appear if the density 
insi de the neutron star e xceeds 2-3 times normal nuclear den- 
sity dGlendenning|1992l) . 

One distinguishes between so-called hybrid stars and the 
strange or selfbound stars. For hybrid stars, quark matter is 
present in a mixe d phase with hadrons and can form a pure 
quark matter core dSchertler et al.ll2000l : lGlendenning & Kettneri 
|2000|) . In the case of a strong first order phase tran- 
sition between the hadronic and the quark phase, hy- 



Bani k & Bandvopadhvavl2003l) . Besides white dwarfs and neu- 
tron stars they constitute a third family of compact stars with 
radii smaller than the ones of neutron stars. If strange quark 
matter is more stable than ordinary nuclear matter so-called 
strange or selfbound stars can form. They are entirely com- 
posed of absolutely stable strange quark matter covered with 
just a thin nucleo n or strangelet crust ( dAlcock et alj 119861 : 
lHaensel et al 111986b. for an introductio n to strange stars see e.g. 
dSchafrner-Bielichll2005l:IWeberll2005l) . The first pure quark star 
was calculated by ItohT dl970» followed up today by a large 
sample of approaches f or quark matter (for an overview see 
ISchaffner-Bielichl d2006l) ). The formation of quark matter can 
take place very early after the formation of the neutron star in the 
supernova explosion, during the proto-neutron star cooling stage 
or it can be delayed by timescales of the order of d ays or years 
depen ding on the mass of the metastable star ( see iDrago et akl 
d2004l) and references therein). Despite their compactness quark 
stars can have masses up to tw o solar mas ses and therefore fit 
well with resent o bservations dOzell [2006) as pointed out by 
lAlford et al.l d2007l) . A X-ray transient with a pulsed component 
of the emission having a frequency f = 1 122 Hz was found r e- 
cently iKaaret et alj ( 120071) . According to lDrago etaT] d2007bl) a 
compact star rotating with that frequency must not only contain 
strange matter but also be a strange star or a quark hybrid star 
due to r-mode instability arguments. 

For arbitrarily high densities and temperatures quarks can 
be treated as a free gas of fermions. However, for neutron 
star interiors qu ark interaction s due to gluon exchange be- 
come important dlwamotolll98ll) . The consequence is that sim- 
ilar to electrons in superconductors, quarks can form Cooper 
pairs with net colour charge which is therefore referred to 
as colour superconductivity. The idea of colour superconduct- 
ing quark mat t er was mentioned alr eady around 30 years ago 
dBairoisl[l977t iBailin & Love1[T98ll) was then revived simul- 
taneously by the work of Rapp et all dl998l) and lAlford et al.l 
dl998[) ( for overvi e w articles abou t colour supe r condu c tivity see 
e.g. iHuand J2005I): lAlford! d200lh: IShovkovvl d2004t): iRischkd 
d2004l) : lRaiagopal & Wilczekld2000l) ). Depending on the strange 
quark mass and the diquark coupling strengths different quark 
pairing patterns are preferred in d ependence of the qu ark chemi- 
cal potential and the temperature. iRtister et al.l d2005l) calculated 
the phase diagram for neutral quark matter with a selfconsistent 
treatment of quark masses. Considering a system of up, down 
and strange quarks where only the u and d quarks are treated 
as massless particles the pairing will only occur between these 
two flavours. The strange quark will not take part in the Cooper 
pairing due to its large mass and the formed colour supercon- 
ducting phase is called the 2SC phase. For quark chemical po- 
tentials which are much larger than the strange quark mass one 
can treat all three flavours as massless. For the highest densities 
with quark chemical potential larger than 400 MeV the up, down 
and strange quarks with all colour charges will participate in the 
Cooper pairing and o ne refers to this cas e as the colour-flav our 
locked (CFL) phase dAlford et alj [19991: iR^petail l200(j) . A 
feature of the CFL phase is the absence of electrons since the 
charge neutrality is provided by the negatively charged strange 
quarks (Raia gopal & Wilczekl200lT) . However for non-zero tem- 
peratures electrons can appear in the CFL phase due to thermal 
electron-positron pair production. For temperatures larger than 
T ~ 10 MeV the metallic CFL phase (m CFL) appears, whe re the 
electron chemical potential is non-zero dRiister et al.l l2004). The 
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so-called gaples s CFL phase (gCFL ) has also a non-zero den- 
sity of electrons dAlford et al.ll2004h . In this phase the pairing of 
strange quarks is suppressed, consequently the pairing gaps for 
the ds-diquark and the us-diquark condensate, respectively, van- 
ish for specific momenta . The gCFL phase has been found to 
have large effects on the cooling behaviour of old neutron stars 
(age < 10 7 years) keeping their temperature warm for a longer 
time period compar ing to neutron stars without the gCFL phase 
dAlford et al.ll2005l) . 

Due to the required quark chemical p otentials (/u a > 300 MeV) 
and low temperatures (T < 100 MeV) dRiister et al.l2005l) colour 
superconductivity is assumed to be realized in the interior of 
strange and hybrid stars. Also from the observational point of 
view there are arguments for it s exis t ence. The high pulsar mass 
of 2.10 +O.28M measured bv lOzell (120061) indicates a very stiff 
equation of state which would result in too fast neutron star cool- 
ing. This fact would favour hybrid stars with colo ur supercon- 
ductin g quark matter cores as was pointed out by iKlahn et"aT] 
(l2006h . Long gamma ray bursts with long quiescent times can be 
explained by a transition from hadronic ma tter to a 2SC phas e 
and then from the 2SC to the CFL phase dDrago et alj|2006l) . 
ISandin & Blaschkel d2007l) study the effect of neutrino trapping 
in new-born quark stars. They find that the cores of newborn 
proto-neutron stars are in the 2SC state and that stable quark star 
solutions with CFL cores exist at low temperatures and neutrino 
chemical potentials. A phase transition from a hadronic phase to 
a quark matter phase or from ungapped to g apped quark matter 
provides a huge amount of energy e oc A 2 /i 2 dAlford et alj|200lt) 
which, when released asy mmetrically in e.g. neutrinos, could ac- 
celera te the neut ron star dDrago et alj|2007 a; Sagert & Pagliaral 
120071) . Recently lNoronha& Shovkovvl d2007l) studied the CFL 
quark phase in the p resence of a strong magnetic field with 
eB/p q 2 % 1 (see also iFukushima & Warringal d2007l) ). The au- 
thors find a smooth cross-over from the CFL to the magnetic 
CFL phase with oscillating pairing gaps and magnetisation for 
increa sing magnetic fields. Furthermore iNoronha & Shovkovvl 
d2007l) argue that the oscillating behaviour of the magnetisation 
due to the Haas-van- Alphen effe ct may cause a large energy re- 
lease during the stellar evolution. ISchmitt et alj (120051) discussed 
the asymmetry in neutrino emission by a stellar core containing 
spin 1 colour superconducting quark matter. Neutrinos are emit- 
ted anisotropically in space due to an asymmetric gap functions. 
Unfortunately, the temperature and therefore the neutrino energy 
for the required phas e to occur is very low a nd leads to too 
small kick velocities. iBerdermann et al.l d2006l) discussed neu- 
trino emission from strange stars in the CFL phase. The mecha- 
nism is based on a beaming of neutrino emission along magnetic 
vortex lines and parity violation of neutrino producing weak in- 
teraction processes in magnetic fields. However, such vortices 
are typically not present in the CFL phase. 
Studying asymmetric neutrino emission from the direct quark 
Urea process: 

d — > u + e~ + v e (3) 
u + e~ — > n + v e (4) 

in a strong magnetic field leads also to a new problem which 
is not present for nucleon Urea processes. For quark mat- 
ter the axial and vector co upling constants are equal and 
simply applying the result of iHorowitz & Piekarewicz] d 19981) 
to quark matter would give a vanishing neutrino asymme- 
try. A s imilar conclus i on ca n be derived from the calcula- 
tions by iDuan & Oianl d2005l) . However, the authors also find 
that an asymmetry is present in the differential reaction rates 
as soon as positrons are involved. Though there has been 



a lot of work done on nucleon and quark direct U r ea pro - 
cesses in strong ma g netic fields ( see iDorofeev et all d!985 ): 
I Duan & Oianl d2005l): lArras & Laild 19991): lLai & Oianl d!998 
bovalldl999T):lBaiko & Yakovlevldl999h : lRiquelme et al.ld2005l 



and llwamotol (1981 



Uaikumar et al.l 



d2006bl) : IChakrabartvl dl996l) : iBandvopadhvav et al.l d 19981) ) 



there are no calculations considering the angular dependence of 
neutrino emissivities. 

In the centre of mass frame the neutrino production from elec- 
tron capture by up quarks is isotropic. Just left handed parti- 
cles are taking part in this process. In the presence of a strong 
magnetic field electrons are forced in the lowest Landau level 
where their spin is polarised opposite to the magnetic field di- 
rection. In the rest frame of the star the quarks have momenta in 
the range of 400 MeV whereas the electron momentum is less 
than 100 MeV. In the boosted frame the momenta of the scat- 
tere d particles are beamed in th e direction of the incident parti- 
cle dBvckling & Kaiantidll973l) . Consequently, the momenta of 
the down quark and the neutrino are beamed in the direction of 
the up quark momentum creating a "neutrino emission cone". 
Positron capture in the centre of mass frame creates a "antineu- 
trino cone" in the same direction. Consequently we expect a po- 
larised neutrino emission along the magnetic axis opposite to 
the field direction. For simplicity we assume that the angle <p be- 
tween the ma gnetic field axis and t he rotational axis is small as 
was found bv lNg & Romanil (|2007|) . 

The paper is organised as follows: In the next section we will de- 
rive analytical and numerical estimates for the degree of electron 
spin polarisations depending on temperature, electron chemical 
potential and magnetic field strength. For a sufficient accelera- 
tion the energy stored in the proto neutron star matter must equal 
the kinetic energy of a pulsar kick. Hence, in Sect. [3] we will 
compute the energy density of quark matter and electrons as well 
as their heat capacities. Subsequently, we will calculate the kick 
velocities analytically and numerically for polarised electrons. 
In Sect. [4] we will discuss neutrino mean free paths in hadronic 
matter and quark matter. Finally, in Sect. [5] we comment on the 
acceleration mechanism for strange stars in the CFL phase and 
close with a summary and an outlook. 

2. Polarisation 

Electrons which are moving in a magnetic field larger than 



eti 



4.4 ■ l0' 1 Gauss 



(5) 



are situated in Landau levels perpendicular to the ma gnetic field 
axis. Their energy is quantised and can be written as dShul'manl 
I2001I) : 

E 2 = m 2 + p? + 2eBT] = m 2 + p 2 + eB[{2v + 1) + 1], (6) 

where the last term in eq. (O corresponds to the kinetic energy 
in the plane. The magnetic field B is pointing in the positive z- 
direction. The Landau Level number 77 is defined by its quantum 
number v and the electron spin s as: 

1 [ +i for n+ 

77 = v + - + s and s = < f (7) 

2 [-j for«_ 

The electron number densities « + and «_ in eq. d7} denote elec- 
trons with spin parallel or anti-parallel to the magnetic field di- 
rection, respectively: 



(2tt) 2 J nmi ,„ Jo 



(8) 
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geB y 
(2tt) 2 ^ 



Jo 



f(E)d P: 



(2n) 2 4^ Jo 



? (^p 2 +ml+2rjeB-ti e )/T + J 



d-Pz- 



(9) 
(10) 



Here g is the degeneracy factor and f(E) is the distribution func- 
tion for fermions. In eq. <j9j the integration over drj is replaced 
by a summation over rj for a small number of Landau Levels. 
For a sufficient high magnetic field all electrons are situated in 
the lowest Landau Level with ?7=0 and spin s = -1/2. In the 
following we will discuss the dependence of the electron spin 
polarisation 



X 



(11) 



on temperature, electron chemical potential and magnetic field. 
As we will see in the next subsection only for zero temperature 
we have a final number of Landau levels. For T + the number 
of Landau levels can become very large. However, the contri- 
bution to the polarisation x of each level decreases rapidly for 
large values of r\ due to the small electron number density and 
can therefore be neglected. FigureQ]shows the numerical results 
for the polarisation of the electron spin for different values of the 
temperature T and the electron chemical potential fi e in depen- 
dence of the magnetic field strengths B. Figure |l(a)| shows the 
polarisation depending on T and p e , where both are scaled by 
the magnetic field. It is assumed that the energy range is much 
larger than the electron mass so that it can be neglected. For 
growing temperature and chemical potential the spin polarisa- 
tion decreases. Such a behaviour is expected, since the additional 
energy helps the electrons to overcome the energy gap of y/2eB 
to the next higher Landau level. However, around T/ VeB ~ 0.1 
and nl yfeB ~ 1.5 the polarisation increases slightly with grow- 
ing temperature. This behaviour is due to a decrease in the den- 
sity of states in the first Landau level which can be calculated 
by 



nil = 1) 



where 



2eB r 



dE 



8(E) 



g(E) = { V^" 



m;-2eE 



E > yjml + 2eB 



0, E < yjm 2 + 2eB 



(12) 



(13) 



and is shown in Fig. [2] for /i = 1.65 VeB. The minimum occurs 
at T ~ 0.15 VeB, which it broadens and shifts to higher temper- 
atures with growing chemical potential. The description of an 
electron gas in the first Landau level is analogous to the case of 
a one dimensional Fermion gas with a mass of yjm~ + 2eB. As 
the temperatures of interest are very small in comparison to the 
electron chemical potential we can make a Sommerfeld expan- 
sion of the number density: 

2eB 



leB C 

= D = 7TS dE 8^ 
(2ny J_ m 



2eB 



(2tt) 2 
2eB 



2^T(l -2 l - 2 ")£(2n)T 2 " 



\d 2 "- l g(E)] 


N 




dE 2 "- 1 






eB 







(14) 



(p 2 - m 2 e - 2eB) 



3/2 
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(a) Spin polarisation for temperature and electron chemical poten- 
tial scaled by the magnetic field 




40 50 60 70 80 90 100 
He [Mev] 

(b) Constraint to fully polarise the electron spin for different mag- 
netic fields in dependence of the temperature and the electron chem- 
ical potential 

Fig. 1. Critical temperature and electron chemical potential p e 
for a given magnetic field strength B up to which electrons are 
fully spin polarised. 



2eB 



(2tt) 2 
2eB 



4 . 2p 2 eB + (eB) 2 



_7 4J« 

30 (,,2 _ ,„2 _ 



42 



(ji 2 -m 2 - 2eB)' 
2p A eB + 6(eB) 2 p z + (eB) 



\3\ 



(K 2 



2eB) 



11/2 



(15) 



The result is also plotted in Fig. |2] and confirms the decrease 
in the electron density of states with temperature. Figure [T(b)| 
shows the constraints on the temperature and the electron chem- 
ical potential to fully polarise the electron spin, that is x — 1 
in dependence of the magnetic field. The temperatures are rang- 
ing up to 10 MeV and ji e up to 100 MeV. The required magnetic 
fields are found to be quite large, between 10 16 G and 10 18 G, but 
lie below the critical value of ~ 1.3 ■ 10 18 Gauss for the stability 
of a magnetised neutron star dLai & Shapiro! [l99~it) . From Fig. 
|l(b)| we see that for temperatures larger than 4 MeV the mag- 
netic fields should be in the range of 10 17 Gauss and higher. Such 
large magnetic fields are not seen in observations but predicted 
to be present during the merger of a binary neutron star system 
dPrice & Ro sswog 2006). Furthermore, the large magnetic field 
is just required to be present in the quark phase and can decrease 
towards the surface of the star. 
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Fig. 2. Electron number density in the first Landau level n(r] — 
1) in eq. ( TT3T > solved numerically and analytically with a 
Sommerfeld expansion for fi = 1 .65 Vefi in the first Landau 
level. 



2.1. Polarisation for vanishing temperature 

The assumption of zero temperature is usually adopted for de- 
scribing cold neutron stars or white dwarfs. The number density 
is given by: 



geB fV" ! -"' ! - 2 ' cii 
(2n) 2 Jo 



(16) 



2 2 

The number of Landau levels is limited to J] max = ^-g where 
jj is the electron chemical potential. For E = fi, the electron 
number densities « T become: 



eB 



eB 



2ffi , 



jj 2 - m 2 — 2veB 



2eB 



(2n) 2 



^ sjp 2 — m 2 — 2veB + -^jp 2 — m 2 



r=l 



(17) 



(18) 



The polarisation is consequently: 

ri- — n+ 1 



X = 



(19) 



+ 1 



For the case (fi 2 - m 2 ) » 2eB the number of occupied Landau 
levels is large and the sum in eq. ( fT9l ) can be transferred back to 
an integration over v: 



2veB 
fi- - m- 



2eB 
p 2 — m 2 



3/2 2 2 

2eB 



2fi 2 - m 2 
3 2eB 

which simplyfies the polarisation to: 



(20) 



(21) 



1 



eB 



X = 



, i 2n 2 -m 2 
^3 2cB T1 



This relation is easy to understand since an increase of the chem- 
ical potential transfers energy to the electrons. Hence, they can 
overcome the energy gap of V2eB to a new Landau Level with 
77 > where both spin directions are possible which decreases 
the value for the spin polarisation^. A look at Fig. |3(a)| where 
the analytical and numerical results are plotted, shows a very 
good agreement between these two approaches. The steps which 
appear in the polarisation curve correspond to the opening of the 
next Landau level. For sufficiently high chemical potentials the 
number of Landau levels becomes so high that the effects from 
the energy quantisation become negligible and the steps disap- 
pear. 

2.2. Polarisation for the non degenerate, massless case 

(fx,m <k V2e5« T) 

The situation which will be discused here can occur in hot proto- 
neutron stars. Within appr oximately 15 s after the supernova ex- 
plosion dPons et al.lll999l) highly energetic neutrinos propagate 
from the proto-neutron star interior outwards. Due to their high 
energy and the large densities the neutrino mean free path is very 
small. While passing through the neutron st ar they heat the inte- 
rior to temperatures of more than 50 MeV dLattimer & Prakashl 
|2004|) . Hence, effects from magnetic fields smaller than 10 17 G 
can be ignored for such high temperatures. For magnetic fields 
in the range of 10 16 Gauss the condition of V2eB «c T holds for 
temperatures larger than 10 MeV. The proto-neutron s tar temper- 
ature drops to T <1 MeV very quickly within 50 s dPons et al.l 
1 19991) and the neutron star becomes transparent to neutrinos. The 
electron mass and chemical potential are assumed to be negli- 
gible in the following. As V2eZ? <K T the number of Landau 
Levels will be very high and we can again replace the summa- 
tion over 77 in eq. ( fTOb by an integration: 

1 



(2^)2 P e E/T + l 



1=0 



eB 
(2 



•b r° r 

^Jo dP J 



1 



e E ' T + 1 ' 



(23) 



(24) 



With /i, m <k y/2eB we have: 

n - = 7T~a ( ' d P\ f dj l 1 ■ 

(2?r) J L Jo e Vp 2+2 ';«B/r + 1 . 

which can be solved with x = ^p 2 + 2rjeB/T: 

3<r(3)r 3 

" 8tt 2 ' 

where £(x) is the Zeta-function with f(3) ~1.202. The number 
density n + is: 



(25) 



(2tt) 2 dP e y/p 2 +2 n eB/T + l 

3£(3)T 3 eB 



a* -4^ rln(2) - 

and the polarisation becomes: 



(22) X = 



1 



ln(2)eB 

n- +n + ~ TTim _ 1 ~ 3r 2 f(3) 



as T 2 » 2eB. 



(26) 
(27) 

(28) 
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From eq. (|28l l we see that for an electron gas with fi,m <k 
y/2eB <K T the polarisation x is proportional to eB/T 2 . The 
higher the temperature, the more Landau Levels are accessible 
for the electrons and the electron spin polarisation decreases. A 
large magnetic field forces the electrons into lower Landau lev- 
els. Figure [3(b)| shows again a good agreement between the ana- 
lytical and numerical approach. The values for^ are very small 
compared to the previous case (Fig. |3(a)] >. The number of Landau 
levels is large and the polarisation goes to zero for high temper- 
atures. The large number of Landau levels is also responsible for 
the absence of steps. 

2.3. Polarisation for the non-degenerate, massive case 

(H = 0, m » T » V2elj 

In this case we will again discuss high temperatures and weak 
magnetic fields. The electron mass m is now assumed to be 
larger than the temperature and the magnetic field. A star with 
such conditions could be an old quark star with a colour super- 
conducting quark phase where the electron chemical potential is 
vanishingly small (e.g. quark matter in the CFL phase). For m + 
the electron number densities have the following form: 

(29) 



(2 



sj 'jr+trr+lqeB/T , J 



Applying the Taylor-expansion for small x = 2eB/T 2 gives: 



A 2 



rjeB 



— r + nx =i 1 

T 2 1 T AT 



where 
A= yfi 2 



+ m L 



(30) 



(31) 



The electron number density n+ can be derived by replacing the 
summation for weak fields with an integration over tj: 



eB 
(2 
eB 



drj- 



1 



1 



T^ AT +1 



eB r 

eB r 
(2nfJm 



dp 

> 

dA 

> 

dA 



AT 
~e~B 



logfl 



+ e 



-A/T\ 



e A ' T + 1 
1 



A 2 T 



log[l 



with 



I ' i i '' ,B 

J{) er at +1 



drj- 



eB VA 2 - m 2 
A 1 

AT 



A/r + 1 
-A/rj 



+ e 



eB 



logfl 



+ e 



-A/T\ 



(32) 
(33) 
(34) 
(35) 

(36) 



For small x the Taylor expansion gives log(x + 1) ~ x and there- 
fore for m » T one finds: 

A 2 T A 



eB r° 

ik 2 ( m2K ° 



eBy/A 2 ' 



m 2 VA 2 - m 2 



-AjT 



m' 




m 


\ eBm 


m 




+ mTK x 








T 




J 


)~(2n) 2Kl 


T 



(37) 



(38) 



where K\ and Kq are the modified Bessel functions of the sec- 
ond kind. T he Bessel functions K n (x) hav e the following approx- 
imate form dAbramowitz & Stegunlll984h : 



K n (x) 



I2nx 



(39) 



for x » n. One should notice here that the index n in K n (x) 
does not appear in the approximation anymore. With the number 
density 







m 




+ mTK x 




7\ 




T 



T 

n_ — 

(2n) 

the polarisation^ results in: 
1 



X 



2T mKg+TK] 
eB K\ 



- 1 



which becomes for m » T using equation d39l : 
1 eB eB 



X = 



2mT 



(40) 



(41) 



(42) 



2(m + T) T 

Analytical and numerical results shown in Fig. |3(c)| exhibit a 
small difference from each other. The reason for this deviation is 
the growing temperature so that at a certain point our assumption 
of m » T is not fulfilled any more. Consequently, the exact 
numerical result and the approximated analytical solution start 
to differ from each other. 

2.4. Polarisation for large magnetic fields (fx, m, T «; ^2eB) 

In this section the magnetic field strength is chosen to be much 
larger than the temperature, the chemical potential as well as the 
electron mass. This situation is in principle the case discussed in 
Sect. 12. ll for large magnetic fields and a non- vanishing temper- 
ature. The number density integral which we have to solve here 
is: 



eB 



H * = (2n) 2 ^ Jo 



n 



dp- 



1 



where: 



eB 



00 /-CO 

V f 



(2*) 2 ^J 



dp- 



1 



eB 



eB 



T ln(2) + 



T ln(2) + 



y/p 2 +2?ieB/T 

v r 



(43) 



(44) 



+ 1 



xTdx 



y> r 

~T J \l2neBIT 



.2 _ 2 n eB 



J2 J 



zdx 



, (45) 



with x = yjp 2 + 2r]eB/T and ^2j]eB/T » 1. The integration 
over x gives: 



eB 



eB 



T ln(2) + Y y/2neBKi 



7=1 



\2neB 



(46) 



T ln(2) + 



(27T) 2 

eBT\n{2) 2eB 



XCO 
drj ^IrjeB 



(2tt) 2 (2^) 2 



K 2 



2r]eBK 1 
2e~B 



^2rjeB 



(47) 



where we replaced the summation over rj again with an integral 
over dr]. With 



2eB 



K 2 



l2eB 



(48) 
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the polarisation becomes 
1 

X = 




(49) 



(50) 



Unfortunately it is quite difficult to approximate the polarisation 
in eq. d50b to a closer expression like the ones in eqs. d22b , (f28b or 
d42l) . In Fig. |3(d)| one can see a sizable deviation of the analytical 
form from the numerical result starting at y/2eB ~ 2T and until 
V2eB ~ IT. The reason might be the transition from eq. (|46*T i 
to (|47] i, where we changed the summation to an integration over 
the Landau levels. The integrand in eq. (|47| | decreases very fast 
with yj2ej]B/T, consequently only the first Landau levels are 
significant. The sum over the Landau levels is always larger than 
the integral over tj. As both, the result of the integration and the 
summation, are standing in the denominator of the polarisation 
(see eq. d50i l) the value forx is larger for the integration, as seen 
in Fig. |3(d)| A solution to this problem might be a mathematical 
more correct way for the transition from the su m to the integral, 
i.e. the Euler-MacLaurin formalism as used bv lSuh & Mafhewsl 
d2000t) . In table[T]we give a summary of all the analytic solutions 
for the polarisation discussed here. 

3. Calculation of Pulsar velocities 

3.1. Heat Capacity 

The specific heat capacity per volume c t . can be calculated from 
the internal energy density e by: 



c v 



(51) 



The energy density of non-interacting relativistic fermions is 
given by: 



4 120 5 ; 8tt 2 



(52) 



with g being the degeneracy factor. Quark interactions can be 
taken into account by perturbative QCD calculations to first 
order in the strong coupling constant a s = g 2 /(4n) (see e.g. 
iGlendenningl (12000^ : 



^ " 120" T 



50a, 
ZLr" 



4 + 8tt 2 



Consequently the heat capacity of quarks becomes: 

50a>. s 



<8e} 



2a s 

71 



I"— I- 



(54) 



21tt 



with the degeneracy factor g = 2x3x3 = 18. Pro to neutron stars 
can have temperatures up to 50MeV with quark chemical poten- 
tials in the range of (400 - 500) MeV. The chemical potentials for 
the electrons are much lower. With \i q » ju e the heat capacities 
for quarks c q are much higher than the ones for the electrons c e . 
The total heat capacity is then to a good approximation: 



I 1 ~ ~~ 



(55) 



3.2. Pulsar acceleration 

As already discussed, the kick mechanism under investigation 
is based on anisotropically emitted neutrinos which accelerate 
the neutron star in the opposite direction of their emission. The 
amount of acceleration depends on the polarisation of the elec- 
tron spin and the neutrino momenta, respectively. Assuming a 
negligible neutrino mass we can write for the kick velocity: 



dv : 



_X_ 

M„ 



(56) 



whereas is the neutrino emissivity, R the radius of the neu- 
trino emitting quark phase and x the fraction of polarised neu- 
trinos. Using €qp = -de/dt the velocity of the neutron star will 
depend on the energy density in the following way: 



-ttR 3 — (e(t ) - e(t f )) = -nR 3 —Ae 



km I Ae 

700—* 7 

s \MeVftrT 3 



1 10km/ 



l-4M e 



(57) 
(58) 



In principle the difference in energy densities in eq. d57| i could 
be produced in various ways in quark ma tter due to the ap- 
pearance of colour-supercon ducting phases dDrago et al.ll2007at 
iNoronha & Sho vkovv 2007). In the following we will concen- 
trate on Ae originating from temperature decay. The right-hand 
side in eq. (l56l l can be replaced by an integration over tempera- 
ture using eq. (IBTI i: 



4 
— ; 

3 



--nR 3 —c g dT. 

3 M ns q 



Hence, 



2 r K, M l - 2 -?h lp ''- T ' 2) 

40 ^i_H ! _fiJ±_-fi_L_) 3 

s \400MeV7 VMeV/ UOkm/ 



L4Mg 

Mns 



(59) 

(60) 
(61) 



where Tq and Tj are the initial and the final temperatures, re- 
spectively. In these calculations we neglect any effects from the 
magnetic field on the energy density of quarks and electrons. 
Considering magnetic field strengths in the range of 10 16 - 10 17 
Gauss and temperatures around 1-10 MeV the phase space of 
the electrons will differ from the one of relativistic fermions. 
However, the electron contribution to the total energy density is 
small and quarks can be considered as a free gas of fermions 
due to their large chemical potential. The situation might change 
when considering gapped quark matter. From eq. d59l ) we can 
see that the neutrino emissivity completely drops out of the ex- 
pression for the kick velocity. The crucial input for the acceler- 
ation is the thermal energy stored in the medium. If the emis- 
sivity is small the neutron star cools slower for a longer time 
period. If the emissivity is large, the cooling is quick. The ve- 
locity estimate in eq. doTt shows that temperatures smaller than 
1 MeV will not give observable velocities. A proto-neutron star 
cools down to temperatures be low 1 MeV during the first min- 
utes (lLattimer & Pra kash 2004). Therefore the final temperature 
T f can be neglected. For a s = 0.5, fx q = 400 MeV, M ns = 1.4M G 
and x — 1 > the dependency of the kick velocity on the tem- 
perature and the radius of the neutrino emitting quark phase is 
plotted in Fig. [4] The higher the temperature of the quark phase 
the higher is the neutrino energy, E v = 3.15T, which results in 
smaller radii R for a given velocity. For an initial temperature of 
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H e 2 /(2eB) (2eB) 1/2 /T 
(a) Polarisation of electrons for T = 0. (b) Polarisation for T 2 » 2eB » m,\i. 



analytical result 
numerical result 



j i i i i i i 

25 2 4 6 8 10 12 14 16 18 

(2eB) 1/2 /T 

(d) Polarisation for/i, m <k T <k V2eB. The magnetic field strength B is 
2- 10 14 Gauss and the temperature is T=0. 1 MeV - 2 MeV. The difference 
between the numerical and the analytical results is max. 10%. 

Fig. 3. Comparison of the numerical and analytical results for the electron spin polarisation as discussed in Sects. 12. 111 - 112.41 

Table 1. Summary of analytic solutions for the polarisation as discussed in Sects. 12. 11 - 12. 41 . 



Case 


Polarisation^ 


Approximation 


T = 

T » 2eB » ii, m 
m»T » 2eB,fi e = 

2eB » T » m,fi 


(22.7T Vi-4^ + ij 

(3se-i)" 

(2T mKolm/Tl+TKtlm/T] 

\eB K\ [m/T] l ) 

\ ln(2) A/ 4<-S e j 


3eB/2fi 2 . for fi e » V2e6, in 

<?Sln2/3f(3)r 2 as T 2 » 2eB 
eB/2mT for m » T 




5 10 15 20 

T/(2eB) 1/2 

(c) Polarisation for m » T » V5efi and fi = 0. 



T > 5 M eV a kick of 1000 km s^ 1 is possible for R <10 km. 
From Fig. |l(b)| we find that these temperatures require magnetic 
fields in the range of B ~ 10 G to achieve full polarisation 
for different electron chemical potentials. These magnetic fields 
seem to be very high compared to the observed ones of approx- 
imately B ~ 10 12 G. Nevertheless, it is possible that these high 
magnetic fields are present in the interior of the neutron star and 
might be substantially smaller at the surface. 



4. Neutrino mean free paths 

For hybrid stars or strange stars with a nuclear crust we have 
to consider at least four neutrino interaction processes with 
the medium: The absorption of neutrinos in quark matter (d + 
v e — m + e~) and neutron matter (n + n + v e — m + p + e~) 
as well as the the scattering processes (q + v — >q + v and n 
+ v — >n + v). We will calc ulate these mean free paths follow- 
ing closely llwamotol d 198 lb . Yet, we have to keep in mind that 
Iwamoto did not consider magnetic field effects and assumed de- 
generate electrons which will not be the case for high tempera- 
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Fig. 4. Radius and temperature dependence of the kick velocity 
for fig = 400 MeV, or, = 0.5, M nj = 1.4 M Q and ^ = 1. 



tures. F or high temperatures one can apply the results of lBurrowsl 
d 19801) who arrives at a temperature dependence of the Urea 
emissivity of T 1 instead of the T 6 behaviour for non-degenerate 
electrons. 



4.1. Absorption of non degenerate neutrinos in quark matter 

For the case of absorption of non-degenerate neutrinos in quark 
matter 



d + v e 



u + e 



(62) 



the neutrino mean free path P. can be calculated for a given 
temperature and neutrino energy by the following expression 
dlwamotolll981l) : 

1 4 , , _ lEl + {nT) 2 \ 

-q— = -^a s G F cos 1 & c p F (d)p F (u)p F (e ) I - - T , (63) 

abs \ I 

where G F is the weak coupling constant and the Cabbibo angle 
is cos 2 0=0.948. The Fermi momenta for the electron, the up and 
down quarks p F (i) with i = e~,u,d can be expressed by their 
respective chemical potential using: 



1 



2a s 



in 



p F (u,d) = fl q , 



(64) 



(65) 



P~u,d 

and 

fl e = PF(e) <K flq 

The strong c oupling constan t a s from quark interactions appears 
necessarily dlwamotolll98ll) where as the prefactor 4/7r 4 in eq. 
d63l differs from the one used by llwamotol dl981l) . 16/n 4 , as 
he applied a modified definition of a s = g 2 /l6n. With a s = 
g 2 /4n = 0.5 the neutrino mean free path becomes: 



I", 

abs 



570km (l + e 



Ey/T 



MeV/ \400MeV/ 



(n 2 T 2 + E}r l 



X \ MeV 2 

VMeV/ V 



Hq 



-)" 2 (— )" 2 . 
') VMeV/ 



(66) 
(67) 



UOOMeV/ 

setting E y = 3. 1ST. One should notice that the relation between 
P, and fi e in equation d67] i is not quite intuitive. The absorption 



reaction in eq. d62l ) should be suppressed for a high number of 
electrons present that is for a high electron chemical potential. 
Hence, one would guess that P, becomes smaller with decreas- 
ing p e . In contrast the scaling of r. with temperature is more 
natural. The larger the temperature the higher is the neutrino en- 
ergy which increases the interaction rate. 



4.2. Scattering of non degenerate neutrinos in quark matter 

For the scattering of non degenerate n eutrinos in quark matter 
and neglecting quark-quark interactions llwamotol d 1 98 ll) gives: 



l< 

l SCClt 



20 



1 



C Vi + C Ai n <°~0 



Pf{i) 



20 



(68) 



whereas i is the quark component flavour, is its number density 
and cr = 4m 2 G 2 F /n. The parameters Cyi and Cai are the vector 
and axial coupl ing constants and can be found in the article by 
llwamotol (119811) . The total scattering mean free path is now given 
by: 



U=M,s scat ! 



20 



2C 2 J + 2C 2 ., + C 2 .. + C 2 „ p 2 q o- m e 



'Vd 



'Ad 



'Vu 



1370km(^^r 2 (-^-y 
UOOMeV/ VMeV/ 

40km (4olfcv) 2 (M^v) ' 



(69) 



(70) 



(71) 



using again E v ^ 3.15T. Since the electrons do not play an im- 
portant role for elastic neutrino-quark scattering their chemical 
potential does not occur in the final expression for the mean free 
path. For fixed quark chemical potentials the interaction rate is 
determined by the neutrino energy. 



4.3. Absorption and Scattering of non degenerate neutrinos 
in hadronic matter 

For the absorpti on and scattering in degenerate neutron matter 
we refer again to llwamotol d 19811) : 



r 

scat 



+3gl)n n o- \ E 



E F {n)j) 

i r 1 

—^{\+3g 2 A )p Fn G 2 F Elm n T\ 

look ( Ev \ 2 ( T PFn V 

VMeV/ I MeV/ \330MeV/ 

1Qkm (JLf (^"L-r 1 

VMeV/ \330MeV/ 



(72) 

(73) 

(74) 
(75) 



where the neutron Fermi energy is defined by E F {n) — p 2 Fn /2m„. 
For evaluating the mean free path in eq. (l75l l we set gA — 1 -25 
and again E v = 3.15T. For the neutrino absorption by neutrons 
we use: 



I" 

L abs 



45km 



(y 4 + IQn-y 2 + 9tt 4 ) 



VIOMeV/ Uo 



4 /„. \-2/3 

no) 



(76) 
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Fig. 5. Neutrino mean free paths in quark matter and neutron 
matter for E v = 3.157\ fi q = 400 MeV, a s -0.5, «„ = «o and 
^=10 MeV. 



~ 230km (M^)"(34Sev)" (??) 
with y = E v /T and n B = «„. For E v » T eq. (IToT l gives: 

" 45km (ToSv)"(53SSv)" < 78 > 



4.4. Discussion of the kick velocity considering neutrino 
mean free paths 

At first sight the mean free paths in the equations ( 167) . (1711) . (l75"l l 
and d77b seem to be quite large. In Sect. 13.21 it was found that 
the required temperatures for the observed pulsar kicks should 
be around 5 MeV. For this temperature the mean free paths de- 
crease drastically to the order of ~ 100 m as can be seen from 
tableland Fig. [5] where we plot the neutrino mean free paths 
in dependence of the temperature. More modern calculations of 
the neutrino mean fre e path in colour flav our locked quark mat- 
ter were done e.g. by iReddv et al.l d2003l) . The authors applied 
for the ne utrino mean free p aths in n uclear ma t ter re sults ob- 
tained by IReddv et all dl998l) and by llwamotol dl98ll) for un- 
paired quark matter. The neutrino energy was set to E v -nT and 
the baryon density is «/, = 5«o which corresponds to a quark 
chemical potential of jj ~ 400 MeV. 

For neutrino mea n free paths in unpaired quark matter 
IReddv et al.l (120031) obtain the same results as in table [2] Large 
differences arise comparing the neutrino mean free paths for ab- 
sorption and scattering in neutron matter. For a temperature of 5 
MeV, the authors get for the scattering a mean free path of 200 
m, which is higher than our result (see table 0, and a mean free 
path for neutrino absorption of 2 m. The large difference in the 
neutrino mean free pat h for absorpt i on is d ue to the different pro- 
cesses applied. While IReddv et al.l d 19981) considered the direct 
nucleon Urea process for n# ~ 5«o we applied the modified Urea 
process assuming a neutron crust where the density is not high 
enough for the direct Urea process to occur. For the next discus- 
sions we will combine the neutrino mean free paths for quark 
matter as well as neutron matter using eq. 



1 



total 



(79) 



with i = n (neutron matter) or q (quark matter). The radius and 
temperature dependence of the kick velocity together with the 
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Fig. 6. Neutrino mean free paths in quark matter and neutron 
matter for E v = 3T, jj q = 400 MeV, a s = 0.5, n n = 0.16 ftrT 3 
and n e = 10 MeV with kick velocities for M ns = 1 AM Q and full 
electron spin polarisation. 



total neutrino mean free paths for quark and neutron matter as 
calculated in eq. d79l ) is plotted in Fig. [6] Different neutron star 
velocities for a given initial temperature T and a quark phase 
radius R are shown. The stars are assumed to have a total mass 
of 1.4M . The dependencies of the kick velocity, the neutrino 
mean free paths and the quark phase radius on the quark and 
electron chemical potentials /i ? and /j. e as well as a s are shown 
in Figs. |6]|9] For a temperature of 5 MeV the mean free paths 
in neutron matter and quark matter are in the range of only 100 
m and therefore up to a factor 100 smaller than the quark phase 
radius. For decreasing a s eq. d54l i predicts an increase of the 
heat capacity and consequently of the kick velocities (see eq. 
(loTli). The same holds for the neutrino mean free path for neu- 
trino quark interactions as seen from eqs. (|63l , d64l and (|68l , 
whereas the absorption and scattering of neutrinos in hadronic 
matter remain unchanged (see eq. d75"l l and eq. (fTTli). A value of 
a s ~ 0.3 is normally expected for quark c hemical potentials in 
the range of 1 GeV (lEidelman et al.l 120041) . But even for such 
a low value we see from Fig. [7] that the change of the neutrino 
mean free path in quark matter and the kick velocities are very 
small. In Fig. [8] we use two different quark chemical potentials 
of iu q = 300 MeV and fi q = 500 MeV. From eq. d54} and eq. 
( |6Tb we expect again an increase of the kick velocity but also 
a decrease of neutrino mean free path for fi q = 500 MeV in 
comparison to Fig. [6] where the quark chemical potential was 
set to 400 MeV. For fi q = 300 MeV, the neutrino quark scatter- 
ing and absorption mean free paths are larger but the kick ve- 
locities decrease. Consequently, we applied a combination of a 
small strong coupling constant a s = 0.3 to increase the neutrino 
mean free path and a large quark chemical potential fi q = 500 
MeV to get higher kick velocities. The results are shown in Fig. 
|9(a)| For quark phase radii of approximately 20 km and temper- 
atures of ~0.5 MeV a kick velocity of 100 km s _1 is reachable. 
Of course such a large radius is not realistic for a quark star. 
Consequently neither the variation of the strong coupling con- 
stant a s nor the change of quark chemical potential will give a 
considerable improvement for the final kick velocities. As can 
be seen in Fig. |9(b)| a significant change in the neutrino-quark 
reaction rate can be achieved by lowering the electron chemi- 
cal poten t ial. H owever, the original mean free paths taken from 
llwamotol (119811) were calculated assuming degenerate electrons, 
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Fig. 9. Kick velocities and neutrino mean free paths in quark matter and neutron matter as in Fig. [6] In |9(a)| the strong coupling 
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Table 2. Neutrino mean free paths for absorption and scattering processes in quark matter as well as in neutron matter with «o = 0.16 
frrT 3 and the neutron Fermi momentum pf„ = (3K 2 n„) 330 («„/«o)^ 3 MeV . 



Medium 


Process 


E v = 3.15T,a, = 0.5 


T=5MeV,// c =10MeV 


Quark matter 
Quark matter 
Neutron matter 

Neutron matter 


Absorption 

Scattering 

Absorption 

Scattering 


labs 30km( MeV ) ( 40 0Mev) (mcv) 
Z fl6l ~230kmU v ) ( 3 30Mev) 

l smt ~ 1 km ( Mc v ) ( 340N ;' C v ) 


~ 120 m for yu f ,=400 MeV 
~ 350 m for /i q =400 MeV 

~ 370 m for n„ = no 
~ 130 m for n„ ^ 5«o 
~ 80 m for n„ = no 
~ 50 m for n„ 5«o 



i.e. /j. e »T. Consequently, for T > p. e the electron chemical po- 
tential in eq. (|67] i is replaced by the temperature, as shown by 
iBurrowsl d 19801) for the neutrino emissivity in the direct quark 
Urea process. Electron-positron pair production due to the high 
temperature and their effects on the neutrino interaction rate have 
been ignored. 

The neutrino mean free path in quark matter hits the 100 km 
s _I velocity curve for R ~11 km and T ~ 1.2 MeV. The 200 
km s" 1 line is reached for R ~15 km and T ~ 1 MeV. In these 
cases the corresponding mean free paths in neutron matter are 
in the range of ~ 6 km. Hence, in this simple approximation, 
which has to be verified in a future work, one can construct pul- 
sar kicks with a velocity in the range of 100 km s and possibly 
200 km s . The magnetic field which is required in this case to 
polarise the electrons can be taken from Fig. |l(b)| and is in the 
range of ~ 10 16 Gauss. The condition p e «: T can be realised 
in a colour superconducting quark phase due to the presence of 
strange quarks (the CFL phase does not require any electrons 
for charge neutrality (Raiago pal & Wilczekl200it) ). Note, that in 
unpaired strange quark matter fi e m 2 J4p q ~ 5 - 6 MeV for 
m s = 100 MeV and fx q = 400 - 500 MeV. Strange stars have 
another feature which would boost the maximum possible kick 
velocity. The size of the hadro nic phase for a hybrid s tar can vary 
between 1 and 3 km see e.g. ISchertler et al.l d2000t) . A strange 
star with M ns = 1.5 M can have a nuclear crust of only 100 
m - 500 m thick ness dGlendenning & W eber 1992). Recently, it 
was proposed bv lJaikumar et alj d2006al) that the crust of strange 
stars can be composed of strange quark nuggets and electrons. 
The thickness of such a crust is calculated to be just 50 m. The 
neutrino mean free path is then similar to a bare strange star and 
one can ignore effects from the neutrino interactions, so that the 
final kick is maximal. 



5. Effects of colour-superconductivity 

In the previous sections we saw that it is difficult to evade the 
problem of the high neutrino interaction rates in hot and dense 
matter by varying a s ,/j. e or \x q . Consequently, we have to look 
for another mechanism. The neutrino mean free path is certainly 
effected by strong magnetic fields. In particular electrons are 
bound to Landau levels. The energy separation between two 
Landau levels is V2efi which would be around 30 MeV for 
B = 1 ■ 10 17 Gauss and in the range of 10 MeV for B = 10 16 
Gauss. The quantised energy levels act as a gap and suppress the 
neutrino quark interaction at least in the plane perpendicular to 
B. Another suppression mechanism arises in quark matter if it 
is in a colour superconducting phase. If the temperature is lower 
than a critical temperature pairing of quarks increases the neu- 
trino mean free paths. To suppress neutrino quark interactions in 
general all quarks have to be paired which is the case in the CFL 



phase. Here charge neutrality is provided by the equal amounts 
of down, up and strange quarks which means that the presence of 
electrons is not required and therefore the CFL phase has ji e = 
dRajagopal & Wilczekll200ll) . On the other hand, electrons can 
be present in the system via electron-positron pair production 
at finite temperature. In addition, a metallic CFL (mCFL) phase 
with nonzero fi e can a ppear for /^-stable m atter for temperatures 
larger than ~10 MeV dRiister et alj|2004l) . L arge electron chem - 
ical potentials are found for fi q ~ 400 MeV dRiister et al.ll2005l) . 
Consequently, it would be interesting to study the neutrino emis- 
sion from quark matter in the mCFL phase as the latter provides 
the quark core with high [i e as well as suppressed quark inter- 
actions. Unfortunately, the pairing between quarks decreases the 
quark heat capacity by: 



c? = 9 (i-^L ? 2 r e -w, (80) 

where A is the gap energy dBlaschkeetal.1120011) . For large val- 
ues of A, c q is lowered drastically. Consequently, the electron 
heat capacity 



He T 



2<r3 



+ — STtT 

30 5 



(81) 



ignored in section 13.11 becomes important. Implementing both 
heat capacities in eq. d59l ) gives the following result for the kick 
velocity: 



-7TR 



3 X 



ln-T 2 
60 



-ttR 3 — r 2 9|i 



2a, 



(82) 



2 

+ —j... 
3 M ns 

As was shown bv lRtister etail d2005l) and lBlaschke et al.1 d2005l) 
large values for the gap A as well as for the electron chemical 
potential p. e can arise in the colour superconducting quark core. 
Since we just want to study the behaviour of the neutrino mean 
free paths and the kick size in principal, we will use the fol- 
lowing estimat es: we choos e the q uark chemical poten t ial to b e 
450 MeV (see iRuster et al.l d2005l) and iBlaschke et ail d2005l) ). 
fi e = 100 MeV and A = 100 MeV. The kick velocity and the 
mean free paths are plotted in Fig. [TO] for a s = 0.5, fi q = 450 
MeV with temperatures up to 50 MeV. These values for T are of 
course quite hi gh, but well in the range of p roto-neutron star evo- 
lution (see e.g. lLattimer & Prakashl d2004l) ). Comparing Figs. [6] 
andfTOlwe see the consequence of the decreased quark heat ca- 
pacity. For radii R < 10 km velocities of 1000 km s~' can be 
reached for T > 30 MeV. For increasing temperatures the sup- 
pression factor from the gap, e~ A/r , diminishes and the contribu- 
tion from quark matter becomes dominant again. As all neutrino- 
quark interactions are suppressed the neutrino mean free path is 
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greatly enhanced for T <K A. However, the kick velocity is re- 
duced at high temperatures so that just values for v < 100 km 
s _1 can be reached for free streaming neutrinos and a realistic 
quark phase radii. As can be seen from Fig. [10] a quark phase 
radius of ~14 km is required to accelerate the star to 100 km 
s . Regarding the high temperatures as well as the high elec- 
tron chemical potential discussed here the magnetic fields for 
fully polarised electrons are expected to be extremely high. For 
T > 10 MeV and electron chemical potentials of /i e > 80 MeV 
the magnetic fields must be larger than 10 18 G and therefore 
higher than the critica l value of ~ 1 . 3 ■ 10 18 Gauss for stability 
dLai & Shapirolll99ll) . iReddv et al.1 d2003l) calculated neutrino 
mean free paths in the colour flavour locked phase to be larger 
than 10 km for a temperature of T = 5 MeV. However, they 
found that the interaction with massless bosons which appear 
due to the breaking of baryon number conservation in the CFL 
phase leads to a mean free path which is similar or even shorter 
than the one for normal quark matter. For temperatures of around 
30 MeV this means a very small value for the mean free path in 
the range of cm which would make the acceleration mechanism 
even in gapped quark matter infeasible. From eq. (|57| | we know 
that a large kick velocity requires a sufficient energy release. Up 
till now the neutrino energy was given by the temperature. In the 
case of a phase transition from an unpaired to a gapped quark 
matter phase with the gap A an energy of 3A 2 // 2 /7r 2 would be re- 
leased ( lAlford et a U2001h . If we assumed a neutrino emission of 
one neutrino per quark, the energy per neutrino would be A 2 //u q . 
For the case A 2 //u q ~ 1 MeV the gap would be in the range of 
20 MeV for fi q ~ 450 MeV. Consequently, for low temperatures 
the neutrino mean free path would be large as seen from eq. d68l 
and eq. j63l . At the same time the energy release 3A 2 fi 2 /7T 2 ~ 3.2 
MeV/fm 3 would be sufficient to accelerate the quark star to 1000 
km s _1 if the quark phase radius was in the range of 8 km (see eq. 
d58l)). In this connectio n it is also interesting to study a scenario 
which was proposed bv lDrago et al.l d2007al) . where the authors 
discuss the transformation of a hadronic star into a quark star 
or a hybrid star. They find the formation of a convective layer 
in the case of conversion from ungapped (e.g. 2SC) to gapped 
quark (CFL) matter. An asymmetry in the fast transport of hot 
CFL quark matter to the surface of the hybrid or quark star could 
create large kicks. 



6. Summary and Outlook 

In this work we studied an acceleration mechanism for pulsars 
based on asymmetric neutrino emission from the direct quark 
Urea process in the interior of magnetised proto neutron stars. 
The anisotropy arises due to a strong magnetic field which forces 
the electrons into the lowest Landau level where their spin is 
polarised opposite to the magnetic field direction. To fully spin 
polarise the electrons for temperatures of 1-10 MeV we find the 
required magnetic field strength to be of the order 10 16 - 10 18 
Gauss, depending on the electron chemical potential. The pulsar 
kick velocity depends on the released energy in the quark phase 
and stems in our case from the heat reservoir. For fully spin po- 
larised electrons we find that for an initial temperature of T > 5 
MeV a kick velocity of 1000 km s is possible for a radius of 
R <10 km and a mass of 1.4 M suitable for a strange star. The 
required magnetic fields have to be in the range of B ~ 10 17 G 
to achieve full polarisation for elec tron chemical potenti als of 
fi e < 30 MeV. However, as shown bv lKusenko et al.l (Q~998) high 
neutrino quark interaction rates will wash out the beamed flux 
of neutrinos necessary for the acceleration mechanism. For typ- 
ical properties of quark matter in the interior of a quark star the 
neutrino mean free paths are in the range of 100 m to 800 m 
dlwamotoll 19811) which is, at least for the quark matter core, too 
small. The best results were found for a s = 0.5, /u q = 400 MeV 
and fi e ~ T giving a velocity of 100 km s _1 for quark phase 
radii of R ~ 11 km and a temperature in the range of T ~ 1 
MeV. A kick velocity of 200 km s was reached for R ~ 15 
km and T ~ 1.2 MeV. For both cases the required magnetic 
field is in the range of ~10 16 Gauss. To suppress the neutrino 
quark interaction and to accommodate a small electron chemi- 
cal potentials we investigated in a first study effects from colour 
superconducting quark matter in the CFL phase. The neutrino 
mean free path is enlarged by a factor exp(-A/T) where A is the 
pairing energy of the quarks. Unfortunately, the quark heat ca- 
pacity also decreases exponentially with A/T so that the largest 
reachable kick velocity is again 100 km s _1 for a gap of A = 100 
MeV, n e = 100 MeV, a quark chemical potential ju ? = 450 MeV 
and a quark phase radius of R =14 km. In summary, the mecha- 
nism discussed here can produce neutron star velocities of 1000 
km s _1 very easily if neutrino quark interactions are ignorable. 
However, the small neutrino mean free paths seems to make it 
impossible to reach velocities higher than 100 km s _1 . For the 
neutrino mean free path to be large small temperatures are re- 
quired. At the same time the energy release should be large to 
accelerate the compact star to high velocities. A solution might 
be a phase transi tion from unpaired to gapped quark matter as 
discussed by e.g. lDrago et al.l d2007al) . where the energy release 
is about 3A 2 jU 2 /^ 2 . The energy carried off by neutrinos would 
be ~ A 2 /fi q per neutrino which is in the range of 1 MeV for 
/j q ~ 450 MeV and A ~ 20 MeV and therefore very small. The 
total released energy is very large and can accelerate the compact 
star to 1000 km s _1 if the radius of the converted quark phase 
is around 8 km. In this c onnection it is intere sting to study the 
phenomenon discussed bv lDrago et al.l d2007al) . The authors find 
the creation of a convective layer at the conversion zone which 
could transport the neutrino emitting gapped quark matter to the 
surface of the star and lead to a pulsar kick. In a future work 
we plan to add effects from quark pairing in the CFL phase as 
well as magnetic field effects for the heat capacities and particle 
energy densities. The energy release in neutrinos and therefore 
the kick velocities as well as the neutrino mean free path in the 
neutron star's interior will be studied in more detail. 
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